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We apply the mutual Chern-Sirnons effective theory (Phys. Rev. B 71, 235102) of the doped Mott 
insulator to the study of the so-called spontaneous vortex phase in the low-temperature pseudogap 
region, which is characterized by strong unconventional superconducting fluctuations. An effective 
description for the spontaneous vortex phase is derived from the general mutual Chern-Simons La- 
grangian, based on which the physical properties including the diamagnetism, spin paramagnetism, 
magneto-resistance, and the Nernst coefficient, have been quantitatively calculated. The phase 
boundaries of the spontaneous vortex phase which sits between the onset temperature r„ and the 
superconducting transition temperature Tc, are also determined within the same framework. The 
results are consistent with the experimental measurements of the cuprates. 
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I. INTRODUCTION 

Since the 1986 discovery^, the high-Tc cuprate super- 
conductors have attracted strong interest both theoreti- 
cally and experimentally. However, a well-accepted un- 
derstanding is still elusive after two decades' efforts. The 
main difficulty comes from the strongly correlated nature 
of the electronic dynamics and the complex phenomena 
observed in experiments. Although there have been a lot 
of theoretical proposals available in describing some as- 
pects of the experimental observations, the most impor- 
tant challenge to a microscopic theory of the cuprates is 
how to provide a consistent understanding of the global 
and universal features of the whole phase diagram. 

The mutual Chern-Simons gauge theory proposed in 
Ref. ^ aims at facing this challenge. This theory is based 
on the phase string formalism of the t — J model^i^, in 
which the charge and spin degrees of freedom are both 
described by bosonic fractionalized fields — called holon 
and spinon, respectively. The fermionic statistics of the 
electron is taken into account by the mutual topological 
interaction between each spinon and holon, mediated by 
the mutual Chern-Simons gauge fields. The advantages 
of this approach are i) naturally including two ordered 
phases, i.e., antiferromagnetic (AF) and superconducting 
(SC), in the phase diagram; ii) explicitly incorporating 
the strong mutual influence between the charge and spin 
degrees of freedom via the mutual Chern-Simons gauge 
structure. In this theory, the bosonic spinous form the 
singlet pairing at a characteristuc temperature Tg, whose 
maximum ^ J ~ 1,480 K at half-filling while monoton- 
ically decreases with doping. Physically, Tq stands for 
the temperature scale below which the short range AF 
correlations start to grow from the length scale of the 
lattice constant.^ The regime at T < To is called the 
upper pseudo-gap phase (UPP), which is the "matrix" 
of all the lower temperature phases for the underdoped 
system, including the AF and SC ordered states. With 
each spinon playing the role of a 7r-vortex of the holon 
field and vice versa, the superconductivity is described 
by the holon condensation with the spinon-pair confine- 



ment, and the AF ordered state is depicted by the spinon 
condensation with the holon self- localization^. Here the 
global phase diagram is essentially decided by the mutual 
duality of the mutual Chern-Simons theory. 

In this paper, we will apply this theory to a low- 
temperature phase embedded in the UPP, which is de- 
scribed by the holon condensate with unconfined spinous. 
Physically, such a regime corresponds to an unconven- 
tional SC fluctuational region, in which strong SC fluc- 
tuations are driven by the low- lying spin excitations. 
In other words, such a phase is both a spin liquid and 
a vortex liquid, named by spontaneous vortex phase 
(SVP). A semi-classical mathematical description of such 
a region in the language of the generalized Ginzburg- 
Landau equation has been given in Refs. m The 
main goal of the present work is to provide a math- 
ematically self-contained effective theory, based on the 
mutual Chern-Simons gauge theory, which can work be- 
yond the Ginzburg-Landau formalism. Starting from 
this effective theory, the physical quantities including the 
diamagnetism, spin paramagnetism, magneto-resistance, 
and the Nernst effect are studied, which are compared 
with the experimental results. The phase boundary of the 
SVP in the T — H (temperature-magnetic field) plane is 
decided by the superconducting transition temperature 
and characteristic magnetic field, Tc and Hm, and the 
holon condensation onset temperature and critical mag- 
netic field Ty , Hc2 ■ AH these four quantities and their 
doping dependence are calculated based on the leading 
order approximation in this effective theory, which show 
qualitative consistency with the experiments. 

The remainder of this paper is organized as follows. 
In Sec. II, the effective description for the SVP is de- 
rived from the mutual-Chern-Simons gauge field theory 
of the phase string model. Sec. Ill is contributed to the 
calculation of physical quantities at the mean-field ap- 
proximation level of the effective theory, which include 
the diamagnetism, paramagnetism, magneto-resistance 
and Nernst coefficient. Then the calculation of the phase 
boundary beyond the mean-field approximation is given 
in Sec. IV. Finally, the conclusion and discussions are 
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presented in Sec. V. 



II. MUTUAL CHERN-SIMONS THEORY OF 
SPONTANEOUS VORTEX PHASE 

A. Mutual Chern-Simons Effective Theory for a 
doped Mott Insulator 

The mutual Chern-Simons theory^ is a field theory de- 
scription of a doped Mott insulator. In such a theory, a 



mutual topological interaction between spin and charge 
degrees of freedom, due to the phase string effect^ in 
the t — J model, is simply captured by a mutual-Chern- 
Simons term. The total lattice Euclidean Lagrangian is 
given as follows^ 



L = Lh + Ls + Lcs 



Lcs 



/ {ij) \ I J I 

J2 bl [do - ^<7AU^)] - | E {^^^'^^ h]_„ + h.c) + A blh, - N {I - 5)\ , 



(1) 
(2) 

(3) 
(4) 



in which Lh and Ls describe the dynamics of the matter 
fields — bosonic spinless holon, hj, and bosonic neutral 
spinon, bi^, respectively. The chemical potential A and 
/J, terms are included to enforce the total number con- 
straint of spinon and holon, respectively. In the fi term, 
S denotes the doping concentration and N the total num- 
ber of lattice sites, and the last term in Lh is introduced 
to account for the on-site repulsion between the holons, 
which may be regarded as a softened hard-core condition. 

The two matter fields hj and bia- minimally couple to 
two U(l) gauge fields A^ and A'^, respectively, in Lh 
and Ls- Physically, the mutual-Chern-Simons coupling 
Lcs [Eq.0)] entangles two otherwise independent gauge 
fields, A^ and A'\ to realize the topological constraint 
due to the phase string effect. This may be directly un- 
derstood by considering the equations of motion for the 
temporal components Aq and A^: 



dL 



dAlil) 
dL 



(5) 
(6) 



In other words, the holon (spinon) number field Uj (Ri^) 
determines the gauge-field strength of A* {A^), as if each 
matter particle (holon or spinon) is attached to a ficti- 
tious TT flux tube seen by the different species. Note that 
in the lattice version of Lcs [Eq.Q], da = Aq for the 
spatial components with AQ,A^(i) = A'^{i -|- d) — A^{i) 

and AaAp{I) = A^p{I)—Ap{I—a). Here the indices a and 



(3 will be always used to denote the spatial components 
(a, [3 = X, y) in the present formalism, and the lattice 
gauge fields A^j = A%{I) {J = I ~ a) and Af^- = A^(j) 
(i — j -\- a)^ with the indices i, / standing for a square 
lattice site and its dual lattice site, respectivelyS,. 

In this theory, the order parameter A|^ = 

X^o^ (|^i<Te"*'^'^''jfej_cr^ 7^ in [Eq.©] characterizes 
the short-range bosonic RVB pairing onset at an upper 
pseudogap temperature Tq. At T < Tq, the AF corre- 
lations start to develop;^ concomitant with the bosonic 
RVB pairing condensation, as described by L^ ■ This the- 
ory predicts^ two ordered phases in the lower tempera- 
ture regions at T ^ Tq, namely, the AF long range order 
(AFLRO) phase with spinon condensation and holon self- 
localization at low doping, and the SC phase with holon 
condensation and spinon confinement at higher doping. 

In the backdrop of the RVB pairing and holon conden- 
sation, the SC phase coherence is realized when spinous 
are confined in pairs. The opposite case is also allowed, 
in which the spinons are not confined, namely, single 
spinous are present as free neutral objects. This defines 
the regime known as the spontaneous vortex phsae (SVP) 
or lower pseudogap phase (LPP), previously discussed 
based on the phase string model in the Hamiltonian for- 
malism at a generalized mean-field (Ginzburg-Landau) 
leveliSii 

In the present work, we will reformulate the descrip- 
tion of the SVP based on the above mutual Chern-Simons 
gauge theory. Since the gauge fluctuations beyond the 
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mean-field level can be more faithfully incorporated and 
treated in this Lagrangian approach, the present formal- 
ism will be more suitable for describing the transitions 
from the SVP to, e.g., the SC phase, as well as the non- 
linear effects like magnetic field dependence of magneti- 
zation, etc. 



\Jn^e^^'. Noting that the amplitude fluctuation vlj of 
hi is gapped, with the phase fluctuation 6i as the most 
relevant mode, the Lagrangian Lh in Eq.Q can be ap- 
proximately reexpressed, up to a constant, as 



B. Description of Spontaneous Vortex Phase 

The SVP is defined as the holon condensed phase, 
in which the holong field can be decomposed as hj = 



I I I{a—x^y) 



I(a=x,y) 



in which is the external electromagnetic gauge vec- 
tor and n'^ = 6. Note that in obtaining the last line of 
the above expression, the following Villain approximation 
has been used in the partition function 



when no density fluctuations are allowed. Due to such a 



background "magnetic" fleld B'^ 



the spin state as 



^jcos[ej-e 

X E ' 



^ ^^^^^ 

[e,-0i-a-A'^{I)-eAl(I)-2iiN^(I)f ^rjs^ 



at large 7 or low temperature. Consequently the resulting 
effective holon Lagrangian becomes quadratic in A* . 

Since the mutual Chern-Simons term is linear in ^f, 
while the spinon part is independent of A^^ , one can thus 
integrate out ^4** in the total partition function to obtain 
the effective dual Lagrangian L'^^^i as: 



rh 



u 
2^ 



E 



i.a 

2* E - E ^"'^AiA^A'^ (8) 



governed by Lg will be driven^ into a spin liquid phase 
at finite doping concentration n'* > as opposed to an 
AFLRO ground state at half-filling. 

The effective Lagrangian L^^^^i also describes the vor- 
tices of the holon condensate as "charge" ±2 particles 
interacting with the Maxwell gauge field A'^ via the min- 
imal coupling A'^J™. By contrast, the spinous in Lg 
also act as the source of the gauge field A'^ with a charge 
(7 = ±1. This corresponds to the fact that in the original 
holon language, each vortex with J™'' = ±1 has a phase 
winding ±277, while each spinon carries a half- vortex with 
a phase winding ±7r, known as the spinon- vortex.-S Ener- 
getically one expects that the charge-2 vortices cost more 
than the charge-1 spinon-vortices and thus, to leading 
order approximation, the former excitations can be ne- 
glected unless when the temperature is very close to the 
upper boundary of the SVP where the phase winding 
starts to lose rigidity. 

But before dropping the term of A'^^J™ from L'^^^i, 
one has to be careful about the case when a ±27r vortex 
is bound to a =F7r spinon-vortex, resulting a ±7r vortex 
which is energetically the same as the original spinon- 
vortex. Generally, including both Lg and L'^^^i, the total 
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where Jq™ = e^^AcA^^ G Z is the temporal component 
of the current for the vortices of the holon condensate and 
J™ = —e^^doNfs denotes the spatial component, which 
satisfy the conservation equation OqJq"'^ + Aq J™'' = 0. 

The above procedure leading to Eq. JSJ is similar to 
the standard boson- vortex duality transformation^. Here 
the field strengths of the Maxwell gauge field Ajj are rep- 
resented by B^^ii) = e"'3A„Aj^(i) and E!^{i) = df^A'^{i) - 
AaAQ{i) and the Gaussian fluctuations of B''{i) — irn'^ 
and E'^{i) are determined by the Maxwell terms in Eq. 
(IHl). Note that the integration over Aq{I) results in the 

same constraint as in Eq. (Q: B''(i) = tttIj which will 1. By defining and $i as the rising and lowering 
be further constrained to Trn'' in the limit of m 00 operators of the holon vortices, with the commutation 



charge coupled to A'l^{i) is given by 
2 Jg°''(i). Consequently, the total charge can be = — 1 
if there is an up spinon co-existent with an anti- vortex of 
the holon condensate at site i, i.e., = 1 (n^ = 0) and 
J™''{i) = — 1. Similarly a down spinon can carry a (7'°* = 
+1 total charge by binding with an vortex with J™' = 
t 
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relations [^,,J^°'ij)] = and $1$, = 1, 

0, the creation operators of the above-discussed spinon- 

and 



vortex bound states can be written as feL^" 



Then, at each site, there will be 4 states that carry the 
minimal total topological charge (7*°* = ±1: 

bh |0) = M +) , bl |0) = |z i -) 
6^^$, |0) = M -) , 6l^$l |0) - |z i +) (9) 

where the sign ± denotes that of (7*°* . The energy of each 
state above is composed of the superfluid vortex energy 
and the spinon excitation. Since the superfluid energy 
only depends on the total dual charge, = 1, the 

four types of spinon- vortices in Eq. |(5J) are degenerate in 
energy. In other words, the spinon- vortex states defined 



in Eq. ((Ujl, rather than the two-component spinous 6^0- 
and the holon vortex field J,™"^, describe the true low- 
energy spin/vortex excitations in the SVP. Similar con- 
clusion has been also reached previously in the different 
mean-field (Ginzburg-Landau) approach for the SVP.^ 
Based on the above discussion, by introducing 
= (^b^^l,bii^ij to stand for the two new 

spinon-vortex states in Eq. JHJ and neglect independent 
±27r vortices, the low-lying effective theory of the SVP in 
the mutual Chern-Simons gauge theory description can 
be finally written down as 



Lcff — Ls 



^dual 



(10) 



where 



= J2 - ^TA'S{^)] ha + E "^l [^0 + - 7 E (C^l-.e^'^^-- + 6L6]-.e— ^fe + h.c 



{ij)a 



+ A 



fh 

^dual 



U 
2^ 



E(4 



iabia + b\^h^ 



N{l-5) 



(11) 
(12) 



in which the Zeeman coupling between the spinon and 
external magnetic field B"^ is also included, with g as 
the Lande g-factor and the Bohr magneton. Such 
an effective theory describes 4-flavor spinous minimally 
coupled to the Maxwell gauge field. There will be two 
possible phases as the consequences: a confined phase in 
which all spinon-vortices form short-range neutral pairs 
and a deconfined phase in which they are free. Physically, 
the former corresponds to the superconducting phase as 
discussed in Ref.^ and the latter is a spin liquid and at 
the same time a vortex liquid, which corresponds to the 
SVP to be further explored in this work. 



III. MEAN-FIELD APPROXIMATION AND 
PHYSICAL PROPERTIES 

A. Mean-field Approximation 

The physical properties of the SVP is described by the 
effective Lagrangian (jlOf) . in which Lg in Eq. (|ll|l deter- 
mines the spinon degrees of freedom minimally coupled 
to the gauge field A'^, while the latter is governed by i^uai 
in Eq. (|12() . Note that the external electromagnetic field 
A'^ does not directly couple to the spinous except for the 
Zeeman coupling, indicating that the spinous indeed do 
not carry electromagnetic charge. On the other hand. 



and are coupled by a mutual Chern-Simons term 
in Eq. (|12|l . implying that the spinous do carry super- 
fiuid vortices, which will become clear in the following 
mean-field solution. 

-'^duai indicates the following mean-field solution for 

Ah . 
■ 

B''(^) = 7rf^^ E'^ii)^0 (13) 

while Aq couples to the external magnetic field by 
—{N^AqB'^. To the leading order approximation, the 
fiuctuations of A^ and Aq around A'^ and Aq may be 
neglected here if the temperature is not very close to 
Tc or Ty, namely, the phase boundaries. When T ap- 
proaches Tc, for example, the long-range interaction be- 
tween spinous mediated via the gauge field will be- 
come important and the fluctuation of in L'^^^i can 
no longer be omitted^. We shall leave the discussion of 
the phase boundaries to the next section. In the follow- 
ing, we first focus on the physical consequences of the 
SVP at the mean-field level, where the spinon-vortices 
proliferate and the long-range interaction between them 
is well screened. 

The effective Lagrangian (jlOf) is then reduced to a 
mean-field Lagrangian 

L^/ = Lg{A\Al, A) - tN-A^,B^ (14) 
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which can be diagonahzed to 

constant by the standard BogoHubov transformational 



with 



m 



(15) 



(16) 



and 



/ A + Em ,^ ^ — Em . , 

Urn. = , Vm = Sgn(f„i)W — (17) 



2E„ 



2 En 



The wave function Wma-{i) = ^m<T(*)j ^^"^ 
Wmai'i) satisfies the Schrodinger equation 

By minimizing L^g^ with regard to the parameters A, 
= A" , and A!^ , the fohowing self-consistent equations 



are obtained 



1-6 



(18) 



N 

B<=a2 1 



N 



(20) 



the first two of which self-consistently determine the av- 
erage spinon number and the RVB order parameter, and 
the last equation describes the polarization of the total 
vorticity of the spinon- vortices by the external magnetic 
field. Here and below, (f>o = ^ stands for the flux quanta. 
Once A, A* and Aq are decided self-consistently, physi- 
cal properties of the SVP can be calculated straightfor- 
wardly. 



B. Magnetization 

The total magnetization can be decided by 

1 dF 



Aft, 



V dB" 



(21) 



with F the mean-field free energy F 



and V = Na^d 



the average volume per CUO2 layer {d is the inter- 
layer spacing). From Eqs. p4ll and pi|l . the total 
magnetization can be then expressed explicitly as 



-AA'n 



9^J■B 



l^d ' 



Mr, 



(22) 



in which Mdia oc iA^, Mpara — (S^) stand for the or- 
bital diamagnetism from the vorticies and the paramag- 
netism from the Zeeman coupling, respectively, where 
iAQ is decided by the self-consistent equations l|18|l . 
Compared to the conventional vortex liquid theory, the 
origin of both the diamagnetism and paramagnetism here 
is intrinsically related to the same spin degrees of free- 
dom, which is a unique feature of the present mutual 
Chern-Simons theory. 

The magnetic field and temperature dependence of the 
total magnetization at different doping concentrations as 
well as the diamagnetism part A/dia at (5 = 0.125 are 
shown in Fig. ^ based on the above mean-field theory. 
The contour plot of Mjia in the temperature and doping 
space at i?^ = 2 Tesla is presented in Fig. |21 Note that 
the diamagnetism disappears at 5 — > where the density 
of the condensate vanishes, as well as at 5 — > ccrvb — 
0.25 (see Refjl) where the RVB pairing vanishes with 
the proliferation of the unpaired spinous at T = 0. The 
magnitude of diamagnetism we obtained is comparable to 
the experimental observations-'^^ in the weak-field region, 
but over-estimated under strong field B"^ ^ -ffc2, where 
the upper phase boundary of the SVP is reached and the 
mean-field approximation is not applicable. 

The paramagnetic magnetization determines the spin 
paramagnetic susceptibility Xs = Mp/ B'^ whose temper- 
ature dependence at (5 = 0.125 and B^ = 2 Tesla is 
shown in Fig. |21 Xs exhibits a prominent "spin gap" 
behavior with reducing temperature below 100 K sim- 
iliar to the pscudogap behavior in the zero field limit ». 
Furthermore a Curie-type upturn emerges at lower tem- 
perature within the "spin gap" , which counts for the con- 
tribution from the free moments in the vortex cores. Ac- 
cording to Eq. H18|) . the total vorticity density satisfies 
n+ — n_ = B^a? and at low temperature when the 
thermal excitations of spinous are suppressed, one finds 
rt_ ~ 0,rt+ ~ B'^o?/(j){), which means there are on aver- 
age Us — B'^a? unpaired spinous per site. At weak 



magnetic field B'^ ^ 0o/a 



^ 1 and the unpaired 



spinous are very dilute and thus are nearly-independent 
to each other. Consequently, a Curie- type spin suscep- 
tibility Xc = ns/kBT is expected as shown in Fig. O 
Experimentally, the NMR spin-lattice relaxation rate for 
the in-plane oxygen nuclear spins, which is related to 
Xi5 due to the hyperfine coupling constant, does show a 
Curie-type temperature behavior in the vortex cores in 
the superconducting phas o^^i^^ , before it eventually de- 
creases and vanishes at very low temperature which is 
presumably due to the Kondo screening effect of the free 
moments by the nodal quasiparticles. 
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FIG. 1: The dependence of the total magnetization M = Mdia + Mpara at various doping concentrations: (a) S — 0.078, 
(b) 5 — 0.125, (c) S = 0.188; (d) the diamagnetism Afdia at 5 = 0.125. The parameters in the mean- field theory are chosen as 
a = 5.5i, d = 1.1 A, J = 120meV. 




FIG. 2: The contour plot of the diamagnetism Mdia (in units 
of A/m) under a weak field B" = 2T . 



C. Magneto- Resistivity 



FIG. 3: The spin susceptibility Xa = Mpan^/B'^ for doping 
5 = 0.125 and B'^ — 2T. Inset is the detail of low temperature 
behavior. The dotted line stands for the Curie curve Xc = 
Y^, with Us = B^a^ /4>o. 



The charge resistivity is non- vanishing in the SVP, al- 
though the holon condensation is still present. The origin 
of the dissipation is due to the flow of the spinon- vortices, 
whose response to an external electric field can be seen 
by taking the classical equation of motion for in the 
dual Lagrangians, l(TT)) and ((T^ . ^ttt = 0, which leads 



to the spinon- vortex current: 

(23) 
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In a uniform stationary state, the last two terms vanish 
and thus one obtains 



JSV 



eaf3cE' 



(24) 



after recovering the full units. 

In a similar way, the equation of motion for Al^ based 
on the original Lagrangian ^ leads to the following re- 
lation 



- 



(25) 



where = iej^ is the holon charge current which 

decides the strength of " electric" field for the gauge 
field A\^. 

Since the spinon-vortices directly see in the La- 
grangian Hll|l , there is generally a linear response relation 
between E'* and J"": 

= <,E^p (26) 
where cr^^ denotes the spinon-vortex conductivity. By 
combining (|24l) - (|26|l and noting that the holon current is 
equal to the electric current in the mutual Chern-Simons 
theory, one finally obtains the electric resistivity as fol- 
lows 

Pafi = - (0o/c)^ ea7Cr^^e5/5 

If one uses a simple semi-classical approximation by ex- 
pressing the spinon-vortex conductivity as ct*^ ~ ^Sa[3, 
with Hy and rjs the total vortex number and viscosity, re- 
spectively, then the charge resistivity in the SVP is given 
by 

P = Pxx^ — [ — ] (27) 



The resistivity in (|27|l is similar to the flux-flow re- 
sistivity in a Type II superconductor except that n„ in 
general is not simply proportional to the external mag- 
netic field B'^. Namely, in the SVP the spinon-vortices 
can be spontaneously (thermally) generated with Uy ^ 0, 
such that p ^ even at = 0. The resistivity p(i?'^) 
can be expanded as 

p{B^)=p{0)[l+jB'^'^ + o{B'')] (28) 

where the odd power terms of B'^ vanish due to the sym- 
metry p{B'^) = p{—B^). Suppose that the dependence of 
the viscosity r^s on B'^ is negligible, then the quadratic 
coeflticient 7 can be expressed as 

^ p(i?-) - p(0) n.(i?-)-n.(0) 
^ p(0)B^2 - n„(0)5^2 ■ ^ ^> 

The spinon-vortex density Uy is determined in the 
mean-field approximation by 



-E 



m,fT 



(30) 
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FIG. 4: The magento-resistence coefficient 7 vs temperature 
for the magnetic field B'^ which is either perpendicular or 
transverse to the ab plane. 



With Eqs. and (jSOJ), the coefficient 7^ and 7||, with 
the external magnetic field B'^ perpendicular and paral- 
lel to the 2D plane, respectively, can be calculated nu- 
merically as shown in Fig. 01 An important prediction 
of the present theory, as shown by Fig. 21 is that 7|| 
is comparable to 7_l in the SVP. This is a rather un- 
usual case for a vortex-flow-induced resistivity, since nor- 
mally the in-plane vortices are always created by the per- 
pendicular magnetic field in a Type II superconductor, 
where the vortex-flow-induced resistivity only exhibits 
field-dependent magneto-resistivity for the compoenet of 
B"^ which is perpendicular to the plane. Experimentally, 
the c-axis resistivity shows an insulating behavior in the 
pseudo-gap phase until T '--^ Tc at low doping, which 
implies that the interlayer quantum phase coherence is 
not important. Thus a magnetic held parallel to the ab 
plane is not expected to contribute significantly to the 
in-plane resistivity based on a conventional flux-flow pic- 
ture, which would predict ~ 7|| ^7^. 

But in the present theory, vortices are tied to the free 
spinous in the SVP. Since the latter can be created by 
the Zeeman term with the external magnetic field point- 
ing at any direction, the former can thus be created by 
the in-plane field as well, although the total vorticity of 
the 2D orbital supercurrents still satifies the constraint 
(|20|l in which B"^ should be replaced by B^. Although 
the present mean-field result of p may not be expected 
to be quantitatively accurate in view of possible correc- 
tions from the fluctuations, the existence of an anomalous 
transverse magneto-resistivity with 7|| comparable to j± 
remains a peculiar prediction based on the mutual Chern- 
Simons theory, which is qualitatively consistent with the 
experimental results for the underdoped YBCOiMii. 
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D. Nernst Effect and Vortex Entropy 



The Nernst effect refers to 



a transverse electric field 



y induced by applying a temperature gradient 
The Nernst coefficient is defined as 



(31) 



One important mechanism that can lead to a significant 
Nernst signal is the flux flow in the vortex liquid phase, 
the contribution of which can be calculated by using Eq. 
(|24|l . Here the vorticity current driven by a temperature 
gradient — V^T is 



(32) 



which, combined with Eq. (|24|l . decides the Nernst coef- 
ficient 



TV 



The vortex thermo-conductivity formula 1)32(1 provides 
a systematic way to calculate the Nernst coefficient. As 
in the last section, we can introduce a simple drift ap- 
proximation usually used in a vortex liquid phase to ob- 
tain a leading order estimation of the Nernst effect in 
the SVP, and leave the more microscopic calculation to 
future works. 

If the vortices drift with a velocity under the tem- 
perature gradient, the vortex current JJ" can be approxi- 
mated by — (n+ — n-)vx = 'j^'^x, in which the last 
equality comes from vorticity constraint (|20|1 . and the 
magnetic field is understood as applied in the perpendic- 
ular direction in the following. The velocity Vx can be 
decided by the equation —Sci,VxT = rjsVx, in which is 
the transport entropy per vortex, and r]s is the same vor- 
tex viscosity as in the resistivity formula 1)27(1 . Thus we 
have Vx% — B'^s^/cporjs, which leads to the Nernst signal 



cr]s 



(33) 



According to Eqs. and the viscosity rjs can 
be eliminated by the ratio 



(poUy 



(34) 



which has the dimension of entropy and relates the trans- 
port entropy and the density of vortices to the observ- 
ables Cy and p. For the conventional Abrikosov vortex 
liquid, one has B'^ = (fioUy and thus S^f, = s^, as usually 
used in the analysis of experiments^-^-. However, in the 
present SVP, both n-)_ and n_ are non-vanishing such 
that generally Sif, < s^. In the mean-field theory, 
can be estimated by = S/riy, in which S is the en- 
tropy density of the spinous, and the numerical results 
are presented in Fig. Oat different magnetic fields. Note 
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FIG. 5: The ratio S^/ks ~ (j>oey/pkB vs. temperature at 
5 = 0.125, and various magnetic field strengths. 



that such an estimation does not distinguish the trans- 
port entropy and those which do not contribute to the 
transport (such as the vortex configuration entropy), and 
may generally lead to an overestimate of and thus S^. 
Furthermore, 5"^ should be further reduced due to the 
vortex-pinning effect at low temperature in the supercon- 
ducting phase, which is not included in the present mean- 
field theory. Experimentally, the Nernst effect in the SVP 
has been studied systematicallyti^*iSiiLi2iiS In Ref. 0, 
the Nernst signal and resistivity are measured for the 
underdoped Lai.92Sro 08CUO4, which shows Cy ~ 7/iV/K 
and p ~ 2 X lO'^n ■ m at H = lOT, T = 25K. Taking 
the c-axis lattice constant d ~ 13 A for LSCO^*^, we ob- 
tain ~ O.T/cb, which is about one order of magnitude 
smaller than our estimation. Nevertheless, Fig. pro- 
vides an upper bound of the spinon transport entropy, 
which drives the spinous downward the temperature gra- 
dient. 



IV. PHASE BOUNDARIES OF THE 
SPONTANEOUS VORTEX PHASE 

In the last section, the mean-field approximation of the 
mutual Chern-Simons gauge theory is applied to describe 
the bulk properties of the spontaneous vortex phase. The 
fiuctuations beyond the mean-field solution will become 
important when the phase boundary of the SVP is con- 
sidered. 

The lower phase boundary of the SVP at low temper- 
ature corresponds to the superconducting phase transi- 
tion. With reducing temperature, the number of ther- 
mally excited vortices in the SVP decreases and the 
screening of the vortex interaction becomes weakened. 
Eventually at T ^ a vortex-antivortex confining 
transition will take place and the system becomes SC 
phase coherent, which will be discussed firstly below 



9 



based on the mutual Chern-Simons theory. 



A. Superconducting transition 

In order to study the superconducting phase transi- 
tion, let us start with the general effective Lagrangian 
(|f 0|l of the SVP. For simplicity, we shall consider the 
limit u ^ oo which means that the holon density fluctu- 
ations are not important and can be neglected. In this 
limit the phase transition will be solely driven by vortices 
as expected. Under such an approximation, the spatial 
component A'^ cannot fluctuate and is constrained by 



Jjh _ - h 



(35) 



Thus the only important dynamical variable will be the 
temporal component Aq. 

By integrating out bi^ , bia in Lagrangian HlOfl 



D[bia\D[bi„] ■ exp 



dTL^ 



(36) 



one obtains a classical action for A^ in uniform external 



field z-(V X A'^) = B"" and E'^ 







ia / 



in which e = ^ D[bic]D\bia\eyiY> drLg 

Under the standard RPA approximation, can be 
expanded to the quadratic order of Aq based on Eq. I|1I(I : 

Fs[AlB^] -z4(q = 0)iV„+i5]xq(S'=)4(q)^"(-q) 

(38) 

in which the total vortex number — {ny{i)) with 
(«) = {^Lbia - bla^iirj and the susceptibility Xq 

is defined by 



Xq 



(39) 

which is equal to the static spin susceptibility according 
to the spinon- vortex binding. Note that here the dynamic 
part of AQ{iu!n) for aj„ 7^ is omitted. By further using 
the expansion 



Xq - Xo (1 + Aoq^) 
in the long- wavelength limit, one finds 



(40) 



FM^„B^] ^ -^A^}(q = 0)iV,„ + f ( E K(0)' + Ao E (^"^S 



(41) 



r 



However, by noting that the totoal vortex number Ny 
is an integer and is conserved in the system described by 
Lagrangian (|I1|) . the exact free energy Fg should be a 
periodical function of Ao(«) in the following sense 

exp I ~I3F, [A'$ + ^, S<=] I = exp [Ag, B"] } 

(42) 

Consequently, the correct form of Fg satisfying H42() as 
well as (|41ll in the weak limit of A^i^i) should have the 
compact version 

^ z^^Aj;i?'= + ixoAoE(^-^o)' 

i ia 

-^^"E'^o^l/'^o) (43) 

in which the self-consistent equation (|f 8() for Ny is used, 
and Xo , Aq can be calculated based on Eq. ifTT)) . 

Finally, combining Eqs. H37|l and H43|) . and redefin- 
ing (3AQ{i) — (f>{ri), we find that ScS reduces to a sine- 



Gordon action 

ScS - 



(frl y (V0)^ - 2ycos(t> 



with 



K 



y 



Xo(S^) 
2/3 ■ 



XoAo 



(44) 

(45) 
(46) 



Such a sine-Gordon effective action describes the fluctu- 
ations beyond the mean-field theory at low temperature 
where the holon density fluctuations are negligible and 
holon 27r vortices are not important. 

In the following, we discuss the superconducting phase 
transition based on this action at low temperature, where 
the phase boundary of the SVP is characterized by the 
superconducting temperature Tc and melting magnetic 
field H„i. 

First, consider the case at B*^ = 0. If xo is negligible 
(the suppression of spin fluctuations) and y ~ 0, the 
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Kosterlitz-Thouless (KT) phase transition temperature 
is decided by the universal value 



Stt 4 



(47) 



which reproduces the conventional KT physics of the 
holon condensate when the long-wavelength spin fluctu- 
ations are absent. Note that Tkt in this limit is only 
1/4 of the usual XY model with the same stiffness thu^ 
because the spinon vortices are tt- vortices instead of con- 
ventional 27r vortices and the latter are neglected at low 
temperature. 

However, the spin susceptibility xo; although sup- 
pressed at low temperature in the present pseudogap 
phase, can be greatly enhanced at the temperature com- 
parable to the spin gap energy (c/. Fig. O, which can 
then suppress the phase transition temperature from the 
universal value Tkt if the latter is higher than the char- 
acteristic spin gap scale. To the linear order of y and K , 
the equation of critical line on the K — y plane can be 
written as2i 



y = 2(1 - 8ttK) 



which results in 



Xo(Tc) 



8^Ao 



(48) 



(49) 



This is a self-consistent equation for Tc as Xo{T) is 
strongly temperature dependent as shown in Fig. |31 
The numerical solution of the superconducting transition 
temperature is presented in Fig. |51 whose value is com- 
parable with the experimental results. By contrast, the 
corresponding unrenormalized Tkt is shown in the same 
figure for comparison. It is should be noted that in the 
cuprate superconductors, the superconducting transition 
can deviate from a pure 2D KT transition22i2^ due to 
quantum fluctuations {e.g., the dynamic fluctuations of 
Aq') as well as the interlayer coupling which may lead to 
a 3D critical behavior near Tc- Nevertheless, the present 
mutual Chern-Simons theory clearly illustrates that the 
low-lying spin correlations play a crucial role in deter- 
mining the temperature scale of the phase coherence. 

Finally, the number of spinon-vortices can be signif- 
icantly induced by the magnetic field via the Zeeman 
coupling as discussed in Sec. Ill, which in turn leads to 
an enhancement of the low energy spin fluctuations as in- 
dicated by the low-temperature Curie behavior shown in 
Fig. 121 Then it is expected that Tc get quickly reduced 
by such an enhancement of the low-lying spin fluctua- 
tions via the magnetic-field dependence of xo- Such a 
suppression of Tc by magnetic field is naturally included 
in the self-consistent equation H49I) . The critical ("melt- 
ing") magnetic field Hm defined by Tc{Hm.) = marks 
a quantum KT transition, where the vortices melt and 
proliferate at T = due to quantum fluctuations. Hm is 
plotted as a function of 5 in Fig. [7| with Tc{H„i) — O.OIJ 
in the numerical calculation. 
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FIG. 6: Doping dependence of Tc decided by Eq. 1)49^ and 
r„ decided by Eq. (I50II . in which the dashed line stands for 
bare KT transition temperature Tkt in Eq. 1471 . All the 
calculations are done for th ~ 3J. 
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FIG. 7: Doping dependence of the "melting" magnetic field 
Hm under low temperature T = O.OIJ, calculated from self- 
consistent equation 1491 with th = 3 J. To the right, the label 
of in the units of Tesla (0o/ci^ = 7100. 9T) is used. 
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FIG. 8: Doping dependence of Hc2 decided by Eq. 1501 . cal- 
culated for th = SJ at T = 0.003J. 
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B. Upper Phase Boundary 

The upper phase boundary of the SVP is defined by 
the temperature and magnetic field scales Ty and i?c2j 
respectively, at which the holon condensation disappears 
and the effective theory H10() breaks down. As has been 
discussed in the Ginzburg-Landau description in Ref ^, 
such a phase boundary can be estimated by the core 
touching of spinon vortices, which completely destroy the 
short range phase coherence of the holon condensate. For 
completeness, in the following we provide the estimation 
of Tv and Hc2 in the effective theory ()10|l , which is quali- 
tatively consistent with the Ginzburg-Landau result but 
quantitatively different. 

The core touching condition can be simply written as^ 

Uy^d (50) 

which means the number of vortices is comparable to 
that of holons. Since we are now focusing on the upper 
phase boundary where the spinon- vortices proliferate, the 
phase fluctuations discussed in Sec. IV A is no longer 
important. Thus we can use the mean- field equation l|3U|) 
of Sec. Ill C to determine the number of spinon- vortices. 
The temperature scale Ty and magnetic field scale Hc2 
of the upper phase boundary are numerically obtained as 
shown in Fig. |Sland|Sl respectively. The value of Ty is 
reasonable as compared to the experiment, whearas Hc2 
is about one order of magnitude largei^, possibly due to 
the omission of the holon amplitude fiuctuations in the 
dual theory (|10|l . which are crucial near the upper phase 
boundary. 

V. CONCLUSION AND DISCUSSIONS 

In conclusion, we have constructed an effective field 
theory description of the spontaneous vortex phase, 
based on the mutual Chern-Simons theory of the doped 
Mott insulator. By introducing a dual transformation 
in the holon Lagrangian, the spontaneous vortex phase, 
which corresponds to a holon condensate without the su- 
perconducting phase coherence, is described as a spinon- 
vortex liquid, in which strong superconducting fluctua- 
tions are determined by the quantum spinon-vortex dy- 
namics. Consequently, both the residual diamagnetism 
and spin paramagnetism can be calculated from the same 
spin free energy under mean-field approximation, and the 
entropy of each vortex can be estimated in this theory, 
which are compared with the results of Nernst experi- 
ments. The key difference between the SVP described 
in the mutual Chern-Simons theory and the other pro- 
posals of vortex liquid phase is the quantum and spin- 
ful nature of vortices, which leads to a closed relation 
between charge and spin properties. As shown in Sec. 



Ill C, the magneto-resistance under transverse magnetic 
field {H II ab plane) is of comparable size as that under 
the perpendicular field H± , which shows that even the in- 
plane magnetic-field can change the vortex number due 
to the spin Zeeman effect. Another prediction from the 
spinon-vortex picture is the existence of spin Hall effect 
in the SVP, which suggests that the transverse vortex 
flow under external electric field and magnetic field 
Bf carries a spin current, as has been discussed in Ref. 

m 

When considering the physics in the neighborhood of 
superconducting phase transition, fluctuations beyond 
the mean-fleld approximation should be included, which 
leads to a KT type critical theory with coupling con- 
stants decided by the spinon-vortex correlation functions. 
The transition temperature Tc calculated from this crit- 
ical theory is shown to be comparable to the experi- 
ments. The magnetic field i?„i needed to kill supercon- 
ductivity can be calculated in a similar way. The upper 
phase boundary Ty and Hc2 are crossover temperature 
and magnetic field scales where the vortex number be- 
comes comparable to the holon number, and the short 
range phase coherence of holons is finally destroyed. In 
the global phase diagram of mutual Chern-Simons the- 
ory, the SVP is a wide fluctuation region on top of the 
superconductivity dome, both of which are embedded in 
the upper pseudogap phase with short-range antiferro- 
magnetic correlations^ 

There are also several issues that are not included in 
the present theory: i) the amplitude fluctuation of holon 
is not included, so that the suppression of holon super- 
fluidity under strong magnetic fleld is underestimated, 
which leads to an overestimation of high-field diamag- 
netism and Hc2- ii) the fermionic quasiparticle is de- 
scribed as a bound state of holon and spinon, which is 
well-defined at low energy, long wavelength regime when 
spinon excitation is gapped. But in the SVP the spinon 
is unconfincd and the quasiparticle will decay into spinon 
and holon, just like the usual RPA collective mode merg- 
ing into the particle-hole continuum. Consequently, the 
contribution of quasiparticle is not important in the SVP, 
although it is important for some low energy features in 
the superconducting phase. 

In the present work, the onset doping of superconduc- 
tivity is (5 = 0, which is a consequence of ignoring quan- 
tum fluctuations of the holon density. Future work is 
needed to derive a more accurate critical theory for the 
superconducting phase transition — both the quantum 
one in the doping axis and the classical one in the tem- 
perature axis. 
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